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Abstract 

We reelaborate on a general method for diagonalizing a wide class of nonlinear 
Hamiltonians describing different quantum optical models. This method makes use 
of a nonlinear deformation of the usual su(2) algebra and when some physical pa- 
rameter, dictated by the particular model under consideration, becomes small, it 
gives a diagonal effective Hamiltonian that describes correctly the dynamics for ar- 
bitrary states and long times. We extend the technique to A^-level atomic systems 
interacting with quantum fields, finding the corresponding effective Hamiltonians 
when the condition of A;-photon resonance is fulfilled. 

PACS number (s): 03.65.Bz, 03.65. Ga, 42.50.Hz, 42.50. Fx 



1 Introduction 



Many nonlinear quantum optical processes, such as fc-th harmonic generation, 
fc-wave mixing or generalized Dicke models, can be described in terms of oper- 
ators from a nonlinear or deformed su(2) algebra [1-8]. The emergence of this 
structure allows us to handle the problem in a close analogy with the usual 
treatment for an angular momentum: in consequence, we get a decomposition 
of the Hilbert space into direct sums of invariant subspaces and, therefore, 
the dynamical problem generated by the corresponding Hamiltonian can be 
reduced to the diagonalization of a finite-dimensional matrix. In spite of this 
considerable achievement, convenient analytic expressions for the eigenvalues 
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and eigenstates are unknown, although several approximations have been de- 
vised relying on the specific features of the particular problem under study. 



Recently [9], the method of Lie transformations (see, e.g. Ref. 10 and references 
therein) has been applied to the approximate diagonalization of nonlinear op- 
tical Hamiltonians that can be recast in terms of the generators of polynomial 
deformed algebras. In this paper we further pursue this topic by applying 
the method to the relevant case of multilevel atomic systems interacting with 
single-mode quantized fields. Consequently, we find the corresponding effective 
Hamiltonians when the condition of /c-photon resonance is fulfilled. 

The contents of this paper are organized as follows: after a brief outline of the 
general method of small rotations developed in Ref. [9] , we show how it works 
for multilevel systems and, finally, we give explicit examples for the cases of 
two- and three-photon resonances. 



2 Nonlinecir su(2) dynamics and effective Hamiltonians 

Let us assume that the interaction Hamiltonian of the model under consider- 
ation can be written as 

H,^, = A Xs + gV, (1) 

where 

V^X+ + X_. (2) 

Here g is a coupling constant that can be always chosen as real, A is a param- 
eter usually representing the detuning between frequencies of different subsys- 
tems (although it is not necessary), and the operators X± and maintain 
the first commutation relation of the su(2) algebra 

[Xs,X±] = ±X±, (3) 
while the second one is modified in the following way 

[X+,X_]=P(X3), (4) 

where PiX^) is an arbitrary polynomial function of the diagonal operator X3 
with coefficients that may depend on some integrals of motion Nj. 
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Now suppose that for some physical reason (depending on the particular model 
under consideration) the condition 

£ = I « 1 (5) 



is fulfilled. Then, it is clear that the Hamiltonian (1) is almost diagonal in 
the basis of the eigenstates of X3. In fact, a standard perturbation analysis 
immediately shows that the first-order corrections introduced by the nondi- 
agonal part gV to the eigenvalues of X3 vanish and those of second order are 
proportional to g/A <^ 1. To proceed further we apply to (1) the following 
unitary transformation (which, in fact, is a small nonlinear rotation) 



with 



in such a way that 



C/ = exp(£r), (6) 



T^X+-X_, (7) 



H,s = UHi^,Ul (8) 



Using the standard expansion 

1 



e^Be-^ = B + [A,B] + j^[A, [A, B]] + . . . , (9) 



we obtain 



00 A; 

H^s^AXs + gY: ^'77— rTTad^(l^), (10) 

where adj- is the adjoint operator defined as adr(V^) = [T, V]. In other words, 
the effective Hamiltonian takes the form 

H,s ^AX^-eg V^X^) + e^^- [X+V^^{X^) + V^^{X^)XJ\ 

- e^^ {XlV^^Xs) + V'$(X3)X! + 2V[$(X3)V2$(X3 - 1)]} 
+ 0{e'). (11) 
where ^{X^) — X^X_ is a structural function such that 

P{Xs) = -V$(X3) = $(X3) - $(X3 + 1). (12) 
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By keeping terms up to order e we finally get 



(13) 



and the resulting effective Hamiltonian is diagonal in the basis of eigenstates 
of X3. 

Furthermore, the higher-order contributions have always the form X^f{X^) + 
f{X^)X^ + g{X^). This makes the procedure of removing the off-diagonal 
terms somehow trivial at each step, in the sense that it is always obvious the 
appropriate transformation that should be applied. For example, to eliminate 
terms of the form 

e\Xlf{X^)+f{X^)X% (14) 

it suffices with applying the transformation 

exp(£^+^Tfc/A;), (15) 

with Tfe = X'lfiXs) - /(X3)X^, since the first commutator of A X3 with Tk 
cancels the corresponding term in the Hamiltonian. 

With this approach, the evolution (as well as the spectral) problem is com- 
pletely solved. The method has the virtue of generality, since it is valid for 
any model whose Hamiltonian could be written down in terms of the gener- 
ators of an arbitrary polynomial deformation of su(2), such as, for example, 
three- wave mixing, Dicke model or k-th harmonic generation [9] . Because the 
transformed Hamiltonian is an expansion in the small parameter e, it is clear 
that its eigenvalues coincide with those obtained using the standard perturba- 
tion theory. The advantage of our method lies in obtaining diagonal effective 
Hamiltonians in an operator form, which avoids the necessity of calculating 
cumbersome series of matrix elements. 

This technique also provides corrections to the eigenstates of the Hamilto- 
nian (1). Indeed, from (8) it is easy to realize that the eigenstates of the 
interaction Hamiltonian (1) can be approximated as 

|*^) = [/t|m), (16) 

where \m) is an eigenstate of X3 and U is the corresponding small rotation. 
Since the rotation operators and \m) are time independent, the operator U 
can be applied to \m) in the form of an expansion in e. For example, the 
eigenstate \^m) up to order takes the form 
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|*J = |m)-£(X+-X_)|m) 

+ + X-) - mXs) + $(X3 + l)]}|m). (17) 

This representation is especially advantageous if we construct the space of 
states of the model as a representation space of the deformed su(2) algebra 
using the raising operator [1^3]: 

|m)ocX^|0), (18) 

where |0) is the lowest weight vector fulfilling the standard condition ^-|0) = 
0. 



3 Nonlinear su(N) dynamics and effective Hamiltonians 

The method of approximate diagonalization embodied in Eqs. (1), (6), and 
(13) can be applied also to systems with a more complicated algebraic struc- 
ture. 

In what follows we consider Hamiltonians which can be represented in terms 
of the su(N) algebra. This algebraic structure naturally arises when describing 
atomic systems with relevant energy levels. 

The Hamiltonian governing the evolution of a collection of A identical A^- 
level atoms (for definiteness, we consider the case of a cascade configuration, 
such that Ei < Ej for i < j) interacting with a single-mode quantum field of 
frequency ujf has te form (in units H—l) 

N N-1 

H = ujfo^a + ujjS" + 9j {aS'^^^ + a^^S^J^^) , (19) 
j=i j=i 

where S*** (i = 1, . . . , A^) are the population operators of the i-th level of energy 
uJi, and S^/ = S'*-', SjI = {S^^Y {i < j) describe transitions between levels i 
and j. The operators 5'*-' form the algebra u(N) and satisfy the commutation 
relations 

[S'\ S^"^] = 5imS^^ + djkS'"^. (20) 
Due to the integral of motion 

N 

YS"=A, (21) 
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one can introduce inversion operators (1<J<A^ — 1) 

gn+i = i ^^j+ii+i _ ^ij) _ (22) 

In this way the operators (5"^, Si^~^^) form the algebra su(N). Apart from (21), 
the Hamiltonian (19) possesses the following integral of motion 

JV-l 

iV = ata + ^ /..-Sf +\ (23) 

with II j — j{N — j), that represent the total number of excitations in the 
system. Let us introduce detunings between energy levels of the atomic system 
according to 

A,- = ujj -ui- {j - l)a;/, (24) 

with the trivial condition that Ai = 0. Moreover, let us assume that A^ satisfy 
the following resonant condition 

A^ = 0, (25) 



which means that the field is in a {N — l)-photon resonance with the atomic 
system: uj^ — uji — {N — l)ujf. Then, the Hamiltonian (19) can be rewritten 
as 

H = Ho + i/int, (26) 

with 

Ho = ujfN + uoA 

Hi^t^ho + V, (27) 
where 

N 

N-l 

V=J29j + a^S'J^^). (28) 
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One can observe that the operators {1 < i < j < N) 

^kk ^ gkk^ (29) 

form a polynomial deformation of the algebra u(N), since they satisfy the 
following commutation relations 

[x'\ x''^] = 0, [x!l, x!!:] = o, [xi^ x'J^] = o, 
[x'',x2] = ±x2iSk,-s,,), 

[X^, X!l] = P{X'\ X*^^), (30) 

where P(X", X''^) is a polynomial function of the diagonal operators X" and 
the integral of motion TV. 

According to the general scheme (6)-(13) we introduce the following transfor- 
mation 

C/i = exp(ri), (31) 



with 

Af-l 



^1 = E «i («'^+ - a^S'J^^) , (32) 



where the parameters af^ are 



and are assumed to be small numbers, a^^^ <C 1, which means that the atomic 
transitions are far from the one-photon resonance (A^+i — = ^^j+i — — 
ujf 3> Qj). Since 

[TiM\--V. (34) 



all one-photon transitions are eliminated by (32) and the transformed Hamil- 

7(1) 

'eflf 



tonian — UiHi^<Jjl takes the form 
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-^eff — ^0 + ^diag + ^nondiag 

JV-2 7. A'-fc-l 

+ E7rTW E ^^''{o'^'S^^'^' + a)'^'s^^'^\ (35) 
fc=i j=i 

(k) 

where the effective interaction constants ijjj can be obtained from the recur- 
rence relation 

^r''-<^^f'-4'^^^^ ■ (36) 

with the initial term tjjj^^ = gj. It is easy to see that ipj''^ ~ af^ and thus 

The term /idiag contains only diagonal terms in the atomic operators and de- 
pends on the integral of motion N (or, which is equivalent, depends only on 
the photon-number operator a'^a). This operator /i-diag appears naturally rep- 
resented as an expansion on the small parameter aj^^ whose first term has the 
following form 

1 N-l 

^diag = 2 E 9,af\Si=-'\2a^a + 1) + {S^^\ Si'^'U + 0{a'). (37) 

The essential point is that this diagonal operator cannot be removed from 
the effective Hamiltonian (35). On the contrary, the operator /inondiag contains 
only nondiagonal terms that can be eliminated by rotations of the type (15). In 
this respect, let us note that the price we must pay for eliminating one-photon 
transitions is the generation of all possible (in this atomic configuration) k- 
photon transitions {k — 2,...,N — 1). The possibility of removing some of 
the terms in the sum appeared in (35) strongly depends on the resonance 
conditions. In the next section we shall examine this delicate point by means 
of some examples. 



4 Examples 

4-1 Two-photon resonance. 

Let us consider the interaction of a collection of A identical three-level atoms 
{N = 3) in a cascade configuration with a single-mode quantum field of fre- 
quency ujf and suppose that all the atomic transitions are far from the one- 
photon resonance with the field. Then, these transitions can be eliminated 
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using the transformation (32). Taking into account the two-photon resonance 
condition A3 = (i.e.; uj^ — uji = 2uOf) and keeping terms up to order cxf\ 
we obtain from (35) an effective Hamiltonian for two-photon transitions (the 
nondiagonal part, /inondiag, can be ehminated by an appropriate transformation 
and no contributions of order o'p appear) 

H^J^ = ho + /idiag + li^?\a'Sf + a^'s]"), (38) 



where 

^diag = lj:gj4\Si^^\2a^a + 1) + Si^+'U], (39) 



and the interaction constants, from (33) and (36), result to be 



(1) 9i (1) 92 

A2 A2 

/(2) (1) (1) r,9l92 ..„x 

n -(^2 9i - 92 = -2-^- (40) 

Taking into account that S^^ + S*^^ + S^^ = A and imposing the absence of 
initial population in level 2 (note that the term ho disappears in this case), we 
obtain the effective two-photon Hamiltonian including the dynamical Stark 
shift [11-13] 



9i92 
A2 



A/2) 




+ A^a^a. 



(41) 



4-2 Three-photon resonance 



Let us consider the interaction of a collection of A identical four-level atoms 
(^N = 4) in a cascade configuration with a single-mode quantum field of fre- 
quency cuf and suppose that there are no one- and two- photon resonances with 
the field. 

After eliminating one-photon transitions the transformed Hamiltonian (35) 
(written up to the corresponding order) has the form 
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-^eff^ — ^0 + ^diag + ^nondiag 



where 

^diag = lEgjaf^iSi'^'C^a^a + 1) + + O(a'), (43) 



2 . , 



and the essential part of the nondiagonal operator /inondiag is 

/inondiag = J E ^^.^ (-S^+^-Si^+' + 5i^+^5!!+l) + ^(a^). (44) 



2 . . T 



The interaction constants are defined according to (33) and (36): 



(1) 9i (1) 92 (1) 93 

1 A' 2 A.-Ao' ^ A, 



(2) 2 A2 — A3 (2) 2 A3 — A2 

,(3) ^9l9293 f.f--. 

where, in addition, the resonance condition A4 = (i.e., three-photon reso- 
nance CU4 — Ui = 3a;/) has been imposed. 

According to the general scheme, the second term (representing two-photon 
transitions) in (42) can be removed by the unitary transformation 

ui'^^eMTh, (46) 



with 



ri^) «f - a^'sl'^'), (47) 



where now 



^i+2 — ^-i 



, ,(2) _ ^ 



are also assumed to be smaU parameters; i.e., there are no resonant two-photon 
transitions (Aj+2 — Aj = iOj^2 — cuj — 2a;/ 3> i^^?) and thus of' <^ ocf' . 
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The nondiagonal term /inondiag can be also eliminated (in the same approxi- 
mation) with the following transformation 

C/f =exp[Tf], (49) 



where 

^2^ - ^ E P^J{Sl-^^S'i-^^ + 4^'+'5!!+^), (50) 



and 



(1) 

A, , 1 - A. + A, - A, 



(51) 



Taking into account that + = A and imposing the condition 

of the absence of initial population in levels 2 and 3, we finally obtain the 
effective Hamiltonian describing three-photon resonant transitions 



- i[ai')yi(ata - 1/2)5^^ - afg^{a^a + 3/2)5^^] 

- ^[«fVf ^a^a(a^a - \)S^^ - afi)f{a^a + \){a}a + 1)S^% (52) 

One must note that the two last terms are of order 1/A^ and, thus, can be 
omitted in the effective Hamiltonian, because the interaction term generat- 
ing transitions between levels 1 and 4 is of order 1/A^. Then, the effective 
Hamiltonian takes the form 



£(ata - 1/2) (A/2 - Sl') + + 3/2)(yl/2 + Sf) 



, (53) 



where ,5,14= (5^4 _ ^ii)/2. 



As in the two-photon case, the three-photon Hamiltonian (53) contains the 
dynamical Stark shift (with the corresponding spontaneous contribution [11]). 
Nevertheless, the essential difference between these two cases hes in the fact 
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that in the two-photon case the interaction term and the Stark-shift term are 
of the same order of magnitude, while in the three-photon case the interaction 
term is one order of magnitude lesser than the Stark-shift term. This lead to 
essential differences in the evolution of some observables. In the same way, 
/c-photon effective interaction Hamiltonian can be obtained. 
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